Dynamical study of the X(3915) as a molecular D*D* state in a quark model 
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Considering the coupling of color 1®1 and 8 ® 8 structures, we calculate the energy of the newly 
observed X(3915) as S— wave D* D* state in the Bhaduri, Cohler, and Nogami quark model by the 
Gaussian Expansion Method. Due to the color coupling, the bound state of D*D* with J^^ = 0^"*" 
is found, which is well consonant with the experimental data of the X(3915). The bound state of 
B*B* with J^^ = 0"*"^ and 2+"*" are also predicted in this work. 



PACS numbers: 12.39.Jh, 14.40.Lb, 14.40.Nd 



O 

(N 
in . 
P.. 

< 



^ ■ 

Oh' 
I 

Oh' 

CD 



> 

^' 

O 

o 



X 



I. INTRODUCTION 

Very recently. Belle Collaboration reported the newest 
charmonium like state, the X(3915), which is observed 
in 77 — ujj I'll) with a statistical significance of 7 .la [J- 
[ij. It has the mass M = 3915±3±2 MeV and width 
r = 17 ± 10 ±3 MeV. Belle Collaboration determined the 
X(3915) production rate 

r^^(X(3915))S(X(3915) ^ ujj/tl)) = 61± 17±8 eV(l) 



or 



r^-y(X(3915))i3(X(3915) -> ujJ/i}j) = 18 ± 5 ± 2 eV, (2) 

for = 0+ or 2+, respectively. The X(3915) is 
similar to the previously discovered F(3940) by Belle 
and confirmed by BaBar Collaboration in the process 
B ifr(3940), r(3940) ^ wJ/V' HQ- The mass 
and width of the r(3940) obtained by Belle are: M = 
3943 ± ll(stat) ± 13(syst) MeV, T = 87 ± 22(stat) ± 
26(syst) MeV jfij. However, BaBar obtained M ~ 
3914.6ti^:^(stat)'± 2.0(syst) MeV, T = 2,4+f{st'At) ± 
5(syst) JVieV Q. The difference of the mass and width of 
y(3940) might be the attribute to the large data sample 
used by BaBar (350 fb~^ compared to Belle's 253 fb~^), 
which enabled them to use smaller ujj/il) mass bins in 
their analysis ^2;, 3:]. The r(3940) and X(3915) have same 
mass and width, and are both seen in ujj/il). Hence they 
might be same states pointed out in Refs.[l|-la|. 

The X(3915) as a conventional cc charmonium inter- 
pretation is disfavored _l- 3]- The mass of this state is well 
above the threshold of DD or DD* Q, hence open charm 
decay modes would dominate, while the loJ /ip decay rates 
are essentially negligible, which is same as F(3940) 9']. 
However, according to the production rate obtained by 
Belle, Eqs.dl]) and ©, we obtain the T^j/^ ~ O (1 MeV) 



if we assume F^-y ^ O (1^ 2 keV) which is typical for an 
excited charmonium state [l0l - [l3j . F^j/^, ^ O (1 MeV) 
is an order of magnitude higher than typical rates be- 
tween known charmonium states, which imply that there 
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is intricate structure for the X(3915). In our previous 
study [3], the quantum number and mass of xo(2'^-Po) 
is compatible with the X(3915). However, the calcu- 
lated strong decay width is much larger than experimen- 
tal data. Hence, the assignment of X(3915) to the char- 
monium xo (23Po) state is very unlikely. 

In the previous work, Liu Jil, [3 suggested r(3940) 
is probably a molecular D*D* states with J^'~^' ~ 0++ 
or J^*^ — 2'^'^ in a meson-exchange model. Assuming 
the D*D* bound-state structure, Branz [l3| studied the 
strong y(3940) -> ojJ/ip and radiative ^(3940) -f-f de- 
cay widths in a phenomenological Lagrangian approach. 
Their results are roughly compatible with the experi- 
mental data about y(3940). By the QCD sum rules, 
Zhang 'is'l also obtained the mass M = 3.91±0.11 GeV 
for D*D*, which is consistent with y(3940) reported by 
BaBar. Using SU(3) chiral quark model and solving the 
resonating group method equation, Liu and Zhang [l^ 
also studied molecular D*D* states with color singlet 
meson-meson structure, the bound state of D*D* may 
appear if the vector meson exchange between light quark 
and antiquark is considered. 

Although at present the data for the Ar(3915) and 
y(3940) are still preliminary, and whether the Ar(3915) 
is the same as the y(3940) needs to be identified, it is 
worthwhile to discuss its possible assignments, especially 
in view of the great potential of finding new particles. 
We assume the A"(3915) and y(3940) as same state in 
the following study. In this work we study the Ar(1915) 
with molecular D* D* structure by taking into account of 
the coupling of different color structures in the Bhaduri, 
Cohler and Nogami (BCN) quark model, since the color 
singlet tetraquark state {cq)*{cq)* can be consisted of col- 
orless or colored (eg)* and (cq)* two-quark mesons. The 
role of the hidden color in qq — qq {q = u, d) has been 
discussed by Brink and Stancu |20|, a nd then extended 
to study other four-quark systems 2l|-[23j . 

The BCN quark model was proposed by Bhaduri, 
Cohler and Nogami [^J]. In this model, the interaction 
between quark pair is composed of linear color confine- 
ment, color Coulomb and color magnetic interaction from 
one-gluon-exchange. This model has been widely used to 
stu dy qq m eson, qqq baryon [2^ and four-quark system 
[2ll. I26l - l28j range from light quark u, d to heavy quark b 
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with same set of parameters. The molecular state bound 
by a one-gluon exchange interaction in the BCN quark 
model can be depicted in Fig. [1] One-gluon exchange 
interaction exist between pair of quark-quark or quark- 
antiquark in color structure 8 (X" 8, while no interaction 
between meson-meson in color structure 1 (g) 1. Hence 
there is no bound state in FigUJA). However, the bound 
state is possible within hidden color channel (color struc- 
ture 8® 8, FigUlB)). 




where r^j = jr^ — r^l and Tc.m. is the kinetic energy of 
central-mass of whole system, a, A are the SU(2) Pauli 
and the SU(3) Gell-Mann matrices, respectively. The 
A should be replaced by —A* for the antiquark. The 
universal parameters of this model 0, [27|] are listed in 
Table H 



TABLE I: Parameters of the BCN model. 





m„,d (MeV) 


337 


Quark masses 


rris (MeV) 


600 




rric (MeV) 


1870 




mt (MeV) 


5259 


Confinement 


Oc (MeV fm-^) 


176.738 




A (MeV) 


-171.25 


OGE 




0.390209 




ro (fm) 


0.4545 




(B) 

FIG. 1: The molecule bound state by a one-gluon exchange 
interaction in color structures 1 ® 1 (A) and 8 (g) 8 (B) in the 
BCN quark model. The darkened and open circles represent 
quarks and antiquarks, respectively. 

The numerical method, which is able to provide reli- 
able solutions, is very important for calculating the en- 
ergy of few-body systems. In this work, the spectra of 
normal and exotic mesons are obtained by solving the 
Schodinger equation using a variational method: Gaus- 
sian Expression Method (GEM), which is a high precision 
numerical method for few body system. The detail about 
GEM can be seen in Refs. 23, 29]. 

The paper is organized as follows. In the next section 
we show the Hamiltonian and parameters of BCN quark 
model. Section|llI]is devoted to discuss the wave function 
of Ar(3915) and similar state Xb with c quark replaced 
by b quark. In Section IIV| we present and analyze the 
results obtained by our calculation. Finally, the summary 
of the present work is given in the last section. 



II. HAMILTONIAN 



The Hamiltonian of BCN model takes the form. 



^ - E ( + i ) - Tern. + E (^^^ + ^rf) (3) 

]>i=l 



with 



yG ^ ( _ . ^ 1 ^ 



4 V ^ij Tfl^Tflj '^o'^ij 



V^ = Xr- AJ {-a,n, - A). 



III. WAVE FUNCTION 

The total wave function of four-quark system can be 
written as. 



(6) 



with 



where |^), jry)^, Ix)"^, \'^) Lt represent color singlet, 
isospin with /, spin with S and spacial with angular mo- 
mentum Lt wave functions, respectively. 

Due to the X(3915) is observed in the invariant mass 
spectrum of J/ipuj in 77 J/il}uj [J|, the flavor wave 
functions of the possible molecular states is 

A:(3915) = ^ (£)0*^o* _^ D-*D+*) . (7) 
v2 

The similar state composed of b quark is 
1 



Xb = (B°*B°* + B-*B^ 
V2 ^ 



(8) 



The spatial structure of the above states are pictured in 
Fig. [2] The relative coordinates are defined as following, 

r = ri-r2, R = r3-r4, (9) 

miVi + m9r9 m^ir^ + miTi 

P ^ : : , 

mi -f- 1712 + 1714 

and the coordinate of the center of mass is 
4 4 



(10) 



Rem = E mir^/ E m^. 



(11) 



(5) where rrii is the mass of the ith quark. 
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FIG. 2: The relative coordinate of a molecular state. The 
darkened and open circles represent quarks and antiquarks, 
respectively. 



By GEM, the outer products of space and spin is 

2m,2 ] ./aMs] J12M12 V'plip)] J.h ■ (12) 

The total angular momentum is or 2 suggested by Belle 
Collaboration l]. 

Three relative motion wave functions are written as, 



(13) 



n=l 



N„ 



^lm(R) - ^NN^LR'-e-'^-'' >lm(R) (14) 



JV=1 



= E Co^N^pP^e-'-'-P^Yp^ip) (15) 

Gaussian size parameters are taken as geometric progres- 
sion 



1 77 — 1 / 



(16) 



The expression of C,N,uJa in Eqs. - ((T5|) are similar 
to Eq. (HH). 

The physical state must be color singlet, so the color 
wave functions of the possible molecular state can be con- 
structed by the following color structures. 



= |ll2(»l34), |6) = |8l2®8. 



(17) 



Clearly, it includes a color singlet-singlet and a color 
octet-octet state. The latter is called hidden color 
states by analogy to states which appear in the nucleon- 
nucleon problem [30| . 



IV. NUMERICAL RESULTS AND DISCUSSION 

Using GEM and the universal model parameters listed 
in Table U we also calculate the mass of normal meson 
composed of quark-antiquark by solving the Schrodinger 
equation 



(18) 



with Rayleigh-Ritz variational principle. The results of 
meson spectra are converged with the number of gaus- 
sians rimax — 7, and the relative distance of quark- 
antiquark ranges from 0.1 to 2 fm which is determined 
by the convergence properties of the energies and dis- 
cussed in detail in Ref. The calculated results are 
listed in Table HIl which are in good agreement with the 
experimental data. 



TABLE II: The normal meson spectra (MeV) in BCN model. 
The last column takes from PDG compilation and the 
the ground state of bottom Tjt{ls) is observed by BABAR 
Collaboration from the radiative transition T{3S) — >■ ^rjt Isil] 



Meson 


BCN 


Exp. 


TV 


137.5 


139.57±0.00035 


K 


521.4 


493.677±0.016 


p(770) 


779.6 


775.49±Q.34 


A'*(892) 


907 


896.00±0.25 


a;(782) 


779.6 


782.65±0.12 


(^(1020) 


1018.5 


1019.422±0.02 


7?c(ls) 


3040 


2980. 3±1. 2 




3098 


3096.916±0.011 


D'> 


1886.7 


1864.84±0.17 


D* 


2021.3 


2006.97±0.19 


Ds 


1997 


1968.49±0.34 


Dt 


2102.3 


2112.3±0.5 


B± 


5302 


5279.15±0.31 


B° 


5302 


5279.53±0.33 


B' 


5351.5 


5325.1±0.5 


S° 


5373.1 


5366.3±0.6 


Bl 


5414.5 


5412. 8±1. 3 


%(ls) 


9422.2 


9388.9+2 3 (stat) 


T(ls) 


9439.5 


9460.30±0.26 



Due to the threshold of four-quark system is governed 
by two corresponding meson masses, so a good fit of 
meson spectra, with the same parameters used in four- 
quark calculations, must be the most important criterium 
[2I H, m m il. Therefore, the universal model pa- 
rameters are also employed in this work to study the 
possible molecular state of the X(3915) and Xb- 

Using GEM with the numbers of gaussians a = 
12, n — 7, N ~ 7, and the magnitude of z ranges from 
0.1 to 6 fm, and 0.1 to 2 fm for x and y, respectively, in 
Eas. (fT5)l - (fT5)) . we obtain converged results for the four- 
quark systems which also discussed in detail in Ref. [23( . 
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By taking the coupling of color singlet and octet into ac- 
count, the total energy of the possible molecular state 
mentioned above are listed in Table IIIII 



TABLE III: The energy of the possible molecular states. Eth 
and E45 are the threshold and total energy of molecular state, 
respectively. The 'cc' stands for channel coupling of color 
structures 1 ® 1 and 8 ® 8. 



Molecular state 




Eth 




E4, (MeV) 








(MeV) 


1® 1 


8® 8 


cc 




0++ 


4042.6 


4043.6 


3941.9 


3912.1 




2++ 


4042.6 


4043.6 


4200.3 


4043.6 


3*3* 


0++ 


10703 


10703.4 


10388.74 


10298.5 




2++ 


10703 


10703.4 


10609.83 


10535.1 



The D*D* with J^^ = 0++ is bound and in good 
agreement with the mass of the X(3915) and Ref. [18| . 
while the J^*-^ = 2++ with same flavor wave function is 
unbound. The results can be understood as follows. For 
convenience. The matrix elements of color operator are 
given in Table ITVl where the (j > i = 1,2,3,4) 

are the indices of quark pairs. The matrix elements 
of spin operator for interacting quark pairs are 
((1,2), (3,4), (1,3), (2,4), (1,4), (2,3)) = (1, 1, - 
2, -2, -2, -2)and(l, 1, 1, 1, 1, 1) for J^"^ = 0++ 
and 2++, respectively. For the color singlet- singlet chan- 
nel, because of no interaction between D* and D*, the 
calculated energies of D*D* must converged to the the 
threshold, the sum of the mass of D* and D*] the dis- 
tances between two quarks (or antiquarks) residented in 
the different color-singlets are much larger than that of 
quark and antiquark in the same color-singlet. The re- 
sults are demonstrated in Table Hill and IVl For the hidden 
color channel, according to the confinement interaction 
in Eq. ([S]) and color matrix element listed in Table IIVI 
the interactions between quark- antiquark pairs (1,2) and 
(3,4) are repulsive, anticonfinement appears [34]. How- 
ever, the other four pairwise interactions are still attrac- 
tive which govern the confinement of whole four-quark 
system. The distance between quark-antiquark pair (1,2) 
or (3,4) is the competition result of the interaction be- 
tween particle 1,2 (3,4) and that between 1,3, 1,4, 2,3 
and 2,4. The distances are converged to values given in 
Table |V] and the energy of the four-quark system is also 
stable (see Table [ml). The color magnetic interaction 
are also important in the hidden color channel. For the 
state with J^'~^ = the contributions from color mag- 
netic interaction — (A^ • Aj)(iTi • dj) of each particle-pair 
are all attractive. Hence, in hidden color channel the en- 
ergy of D*D* with J^^ = 0++ is lower than the state 
with J^'-^ — 2++, where the contributions from color 
magnetic interaction for particle-pairs (1,2) and (3,4) are 
attractive while that for other four particle-pair are re- 
pulsive. The channel coupling increases the difference 
further, because the cross matrix elements between color 



singlet-singlet and color octet-octet with J^'^ = 0++ is 
two times larger than it with J^'^ = 2++. 

From the Table Hill one can see that the B*B* with 
J^'-^ = 0++ and 2++ are both bound states. The ap- 
pearance of the state B*B* with J^'~^ = 2++ is due to 
the larger mass of h quark than the c quark, the kinetic 
energy of former is lower than the latter, the repulsive 
color magnetic interactions in B* B* is also small (the in- 
teraction is inversely proportional to the square of quark 
mass) in hidden color channel which leads to the energy 
below the threshold of B*B*. 



TABLE IV: Color matrix elements (Cfc|A^ • X'j\^,) with kJ = 
1,2 and j >i = 1,2,3,4. 







(1,2) 


(3,4) 


(1,3) 


(2,4) 


(1,4) 


(2,3) 


{k,l) = 
(fc,0 = 
(fc,0 = 


(1,1) 
(2,2) 

(1,2) 


16 

2^ 
3 




16 

2^ 
3 






4 
3 

/32" 
V 9 




4 
3 

/32" 
V 9 




14 
3 

/32" 
V 9 




14 
3 

/32" 
V 9 



TABLE V: The distances between each quark pair in the con- 
figuration of D* D* in different color structure. Vij stands for 
y^< r'j- > {j > i = 1, 2, 3) in this form, (unit: fm) 





Color structure 


ri2 


r34 


ri3 


^24 


rii 


'■23 




1® 1 


0.65 


0.65 


6.6 


6.6 


6.6 


6.7 


0++ 


8® 8 


0.7 


0.7 


0.69 


0.69 


0.36 


0.91 




Coupling 


0.71 


0.71 


0.69 


0.69 


0.36 


0.91 




1® 1 


0.65 


0.65 


6.63 


6.63 


6.6 


6.6 


2++ 


8® 8 


0.8 


0.8 


0.78 


0.78 


0.39 


1.03 




Couphng 


0.65 


0.65 


6.62 


6.62 


6.59 


6.64 



V. SUMMARY 

In summary, we have studied the possible molecular 
states of X(3915) and similar state Xb with b quark re- 
placing c quark by taking into consideration of the cou- 
pling of different color structures in BCN model. In the 
quark model, there is no net interaction between two 
color singlet mesons. To form the molecular state, the 
hidden color channel must be introduced. By considering 
the hidden color channel and the coupling between color 
singlet channel and the hidden color channel, bounded 
states D*D* and B*B* are obtained. In our calculation. 
The X(3915) can be interpret as D*D* with J^'^ = 0++, 
which is in agreement with the discussion in Ref. 0, [13- 

Since the experimental data about the X(3915) is pre- 
liminary, further confirmation by other experimental col- 
laborations is necessary. The role of the hidden color can 
be tested in multi-quark system if the X(3915) is iden- 
tified as a molecular state D*D* with J^'^ = 0++. At 
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present, we only have a limited understanding for the 
puzzling state X(3915). A systematical study of the role 
of hidden color in different four-quark configurations will 
help us to understand the effects of hidden color config- 
urations. The work is in progress. 
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